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Abstract 

Thermodynamics is derived rigorously for the ID supersymmetric t- J model 
and its SU(if, 1) generalization with inverse-square exchange. The system 
at low temperature is described in terms of spinons, antispinons, holons and 
antiholons obeying fractional statistics. They are all free and make the spin 
susceptibility independent of electron density, and the charge susceptibility 
independent of magnetization. Thermal spin excitations responsible for the 
entropy of the SXJ(K, 1) model are ascribed to free para-fermions of order 
K - 1. 
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Elementary excitations in one- dimensional electrons consist of collective modes of spin 
and charge with velocities different from each other. This is often referred to as the spin- 
charge separation. In the supersymmetric t-J model with the nearest-neighbor hopping, the 
spin velocity depends on the average electron density n per site [[L|]. In the model with long- 
range interactions, on the contrary, the spin velocity is independent of n, and the charge 
velocity is independent of the magnetization m We refer to this independence as the 

strong spin-charge separation. A question of basic importance is then whether the strong 
spin-charge separation persists at finite temperature T. The purpose of this paper is to derive 
thermodynamics of the model microscopically. We show that the low T thermodynamics 
is determined by elementary excitations obeying fractional statistics, but that the strong 
spin-charge separation breaks down as the temperature increases toward a characteristic 
temperature. 

At zero temperature elementary excitations in the long-range t-J model have been in- 
vestigated in || where free semionic spinons, holons and bosonic antiholons are identified. 
The thermodynamics of the spin chain (Haldane-Shastry model) has been derived with the 
use of the empirical supermultiplet rule J7|. The empirical rule was generalized to the t-J 
model in writing down the basic equation for the free energy [Q. The symmetry behind the 
supermultiplet is now identified as Yangian ||. On the other hand, a simple explanation 
for the supermultiplet follows by deriving the family of long-range lattice models from the 
continuum Sutherland model JTT[ in the limit of large coupling constant p,10| . In [10 



the asymptotic Bethe ansatz (ABA) was used to derive the spectrum of the Sutherland 
model. In this paper we also utilize the limiting procedure, but do not rely on any unproved 
assumption. 

The supersymmetric t-J model in one dimension || is represented in a form 

H = y £ l t ij (P ij + l-2ni), (1) 

i<j 

where rii is the electron number operator at site i. We have introduced a graded permutation 
operator = Y^ a ,/3-^i-^j^P where X" 13 changes the state a at site i to @ with a, f3 being 
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either (vacant) or o = ±1 (occupied by either spin written also as |, j). The sign factor Op 
is — 1 if (3 — and is 1 otherwise. The interaction has the long-range form: = tD(xi—Xj)~ 2 
with D(xi — Xj) = (L/tt) sin[7r(xi — Xj)/L] and t > 0. Here L denotes length of the system, 
and Xi a lattice site. 

We use as an auxiliary a variant of the Sutherland model in the continuum one- 
dimensional space. It is given by 

where My is the exchange operator of coordinates of particles i and j || . The spectrum of 
the model obtained by the ABA [[TtJ has been proven to be exact The ground state 
energy Eq is of 0(A 2 ), and the energy E relative to Eq is given in terms of the distribution 
function v a {k) for the component a with momentum k. Here k is an integer multiple of 
2tt/L with the periodic boundary condition. The result is 

s-^E^W + ^glfc-VK*).^ (3) 

where u(k) = Y^a^aik)- For identical particles with internal degrees of freedom, the sym- 
metry of the wave function leads to M^P^ = — 1 within this Hilbert space. We take the 
limit of large A and m, keeping the ratio t = \/(ma 2 ) fixed. Here a = L/N corresponds to 
the lattice constant. In this limit the N particles crystallize at the lattice points, and the 
first term without A in Eq.(|3D becomes negligible. The second term describes the spectrum 
of the lattice model given by Eq.(|l|) plus the lattice vibration (phonon). In the following 
we take the units such that t = a — 1. It can be shown that the phonon frequency uj q is 
given for positive momentum q by u q = q(ir — q/2) = (it 2 — p 2 )/2 = u(p), where we have 
introduced p = tt — q and w(p) for later convenience. 

The quasi-particle energy e(k) for the component a is given by e(k) = 5E/5v a {k) which 
in fact is independent of a. The velocity (or rapidity) p{k) is defined by p(k) = de(k)/dk 
which tends to tt as k goes to infinity. By further differentiating the rapidity we obtain 



in the thermodynamic limit iV — > oo. 

The distribution functions are determined so as to minimize the thermodynamic poten- 
tial. Namely one has v a (k) = f (e(k) — p a ) and Uo(k) = b (e(k) — po) where / and b are fermi 
and bose distribution functions with chemical potentials p a . Then integration of Eq.([|) with 
respect to e gives 

i(7r 2 - p 2 ) = = Tln(l + vo) - T ]T ln(l - v a \ (5) 

where ir 2 in the leftmost side comes from the boundary condition at e — > oo, and where 
obvious arguments of distribution functions are omitted. The results so far given are the 



same as those obtained in |L0j with use of the ABA, but are quoted here as prerequisite to 
the new results to be given below. 

It is convenient to introduce another distribution function p a (p) in the p space by 
p a (p)dp = v a (k)dk, with the sum rule J2aPa(p) = 1 for each p. Excitations in the t-J model 
can be described in terms of po(p) for charge and p s {p) = p-\ (p) — Pi(p) for spin. At T = 0, 
the distribution functions reduce to step functions: po(p) = 9(p c — p) and p s (p) = 9(p — p s ) 
where p c = 7r(l — n) and p s = 7r(l — m). Here n = + is the average number of electrons 
per site, and m = — the magnetization. Note that p c and p s corresponds to velocities 
of charge and spin, respectively. 

Magnetic and charge susceptibilities describe changes of m and n against the changes 
of the magnetic field h and the electron chemical potential (. These are related to p a by 
Pa — Pq = C + ah. The phonons are not affected by the change and can be disregarded. Let 
us consider the low T case where P = exp(—j3h) 1 and M = exp(—/3() <^ P. Then 
the thermodynamics is determined by excitations near p c and p s . We shall first derive the 
charge susceptibility and introduce 

e c {p) = \{p 2 -^) + 2^\{p 2 -p% (6) 

Then po(p) near p = p c is obtained as 

= bje - fi ) = 1 

P ° P 2 + 6(e-/i ) ^exp^e^ + l 



where /(e — p a ) has been approximated by 1 with the condition M<P< 1. The function 
2po(p) is the distribution function for semionic particles with energy e c (p)/2 |fL3"] . This 
particle is called the holon. We shall discuss the statistics in detail later. 
The density n is given by 

^ n= r%^ir decP ^(J^), (8) 



/o 7r 7r j-oo y ae c 

where we have used the delta-function like character of —dpo/de c in extending the range of 
integration. The quantities 

are calculated to be: J = 1, ii = 0, J 2 = 2tt 2 T 2 /3, I 3 = 12^(3)T 3 . Thus one can perform a 
low T expansion of n. With the use of Eq.(^|) we obtain p(e c ) ~ Pc+tc/Pc—t 2 -/ (2Pc)+ e c/ (2Pc)- 
The charge susceptibility Xc(^) is given by 

which is independent of m. This independence is a signature of the strong spin-charge 
separation at low T. The presence of 0(T 3 ) term in Eq.flTOD makes a difference from the 
standard Sommerfeld expansion. Analysis of the p integral shows that the m-dependence at 
low T enters through an exponentially small parameter exp(— T mix /T) where 

T m ix = ^{P 2 S -P 2 C ) = ^ 2 (2 -n-m){n-m). (11) 

We now turn to the spin susceptibility. In deriving p s near p s , we can set po = 0. Then 
from Eq.(|5|) we obtain with a little manipulation 

p,(p) = 1 (12) 

y/Aexp\j3e 8 (p)] + 1 

where e s (p) is given by 

es(p)= 1 -(^-p 2 )-2h = ^(pl-p 2 ). (13) 



Thus the spin excitation with energy e s (p)/2 also obeys the semionic statistics, and is called 
the spinon The magnetization m is given by integration of p s (p) and the differential 
susceptibility is derived as 

f)m 2 2T 2 

— = Xm (m) = — -[1 + -] + 0(T 3 ). (14) 

oh n 2 (l — m) 37r 2 (l — m) 4 

It should be emphasized that x m (™) has precisely the same functional form as Xc {n). That 

Xm{ m ) is independent of n is another signature of the strong spin-charge separation. The 

spin susceptibility Xs is related to the magnetic susceptibility Xm by Xm = 4=Xs- 

If we take the limit of zero magnetic field first, i.e., h/T <C 1, we have p s {p) = 

/3hexp [Pu(p)/2] . Then the susceptibility is given by 

X m = l~ 2 (l + ^T) + 0{T\ (15) 

which has the 0(T) correction in contrast to Eq.([L4]). The difference comes from the p- linear 
spinon spectrum near p = it. For general temperature, Xs — Xm/4 at h = can be derived 
numerically from 

Xs = I dp ^rrv- (16) 

Figure 1 shows the results for various n. It is clearly seen that the zero temperature limit as 
well as the initial slope is independent of n. This again demonstrates the strong spin-charge 
separation. We remark that there is no logarithmic singularity near T = in contrast to 
the Heisenberg model [0. The absence confirms that the supersymmetric t- J model is the 
fixed-point model for one-dimensional electrons PJBUI5"]. As seen in Fig.l, Xs does come to 
depend on n with increasing T. The breakdown of the strong spin-charge separation already 
begins at T substantially lower than T mix given by Eq. QTTJ) . 

Let us turn to the entropy Stj per site of the t- J model which consists of the boson part 
So, the fermion part S a for each spin, and minus of the phonon part S p h- We first compute 
S given by 

S Q = r^i[(6 + l)ln(6 + l)-61n6], (17) 

JO TT V 
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where one may set v = 2 + b at low T since the dominant contribution comes from p ~ p c . 
Eliminating b in favor of po we get by partial integration 

The fermion part can be derived for arbitrary magnetization. We start with the expression 
NS a = — J2k[ u (T m ^o- + (1 — Vo) ln(l — Vo)], and change the integration variable to e(k) — p a . 
At low T the dominant contribution comes from e(k) — \i a ~ 0. Then we get 

& ^£ rf£(£ _^ (£) (_^),- (19) 

with p CT = de(k)/dk at e(/c) = p; CT . In the case of m > 0, the contribution with p^ = tc is 
the same as S p h and cancels each other, while pj = p s = 7r(l — m). The final result is 

5u ^ | + ^— ) = -jrT ( Xc + Xm) , (20) 

3 V 1 — n 1 — m/ 6 

which also corresponds to the specific heat jT. The result describes a two component 
(spin and charge) Tomonaga-Luttinger liquid, and proves previous conjectures @,[J. In the 
opposite limit of high T, we may neglect the p-dependence of p a (p). Then for fixed n a we 
recover the obvious result S t j = — X^TVlnno- — (1 — n) ln(l — n). 

The fractional statistics of excitations is better understood by generalizing the super- 
symmetry to SU(K, 1). If a particle obeys the exclusion statistics characterized by g, the 
distribution function p(p) obeys the equation 

p h (p)+gp(p) = l, (21) 

where p h {p) = wp{p) is the distribution function of the anti-particle (hole) |T5[ . The weight 
factor w is equal to 1/vq for the charge component, and to (l—v a ) / v a for the spin component 
a — 1, . . . , K in the SU(K, 1) model. The particle-hole duality P,[13|j becomes apparent if 
one divides both sides of Eq.([21|) by g 0) and makes a rescaling dp — > dp/ g. Then it is 
seen that the anti-particle follows the 1/g statistics. The energy of the anti-particle is —1/g 
times that of the particle with the same p. 
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The holon and antiholon distributions near p = p c are given by po — bj (K + b) and 
Po(p) = l/(K + b) at low T with M < P < 1. Because of the relation p (p) + Kp^(p) = 1, 
the antiholon obeys the if statistics without the rescaling of dp. The energy of the antiholon 
is given by — e c (p). The holon on the other hand obeys the 1/K statistics and p needs the 
rescaling dp — > dp/K. Thus the compensating factor K = 2 has appeared as 2p (p) for 
holons in the SU(2,1) t- J model. 

Similarly p a near p = p a for a = 1, . . . , if is rewritten as p Q = v a / (a — 1 + u a ) where 
the occupation n\ > n 2 > . . . > tik > has been assumed. Then we get 

a Pa (p) + (a-l)p h a (p) = l (22) 

where p„(p) is the distribution function of the anti-particle (spinon). Dividing Eq.(|22"D by 
a we obtain the statistics of the spinon as (a — l)/a with the rescaling dp — > dp/ a. The 
anti-particle of the spinon (antispinon) correspondingly follows the a/(a — l) statistics with 
the rescaling dp — ► dp /(a — 1). The antispinon is first identified in this paper. It has not 
been noticed in |J probably because the singlet ground state with p a = n for all a has no 
antispinon. In the special case of a = 2 the statistics of spinon is reduced to 1/2, i.e. to the 
semionic one. The rescaling explains why the spinon spectrum || is periodic in 7r instead of 
2tt. 

For thermal excitations a description different from the exclusion statistics can be more 
convenient. As is well known the entropy of the X-Y chain can be understood most easily 
in terms of free fermions introduced by the Jordan- Wigner transformation. For spin com- 
ponents more than 2, ref. [0 has introduced para-fermions for another long-range model 
with harmonic confinement potential. We now explore this type of description in the present 
model. For simplicity we consider the case of n — 1 without magnetic field. Setting 6 = 
in Eq.(|3p we obtain u a — 1 — exp(— /3u/K) with u = (tt 2 — p 2 )/2 for a = 1, . . . , K. Then 
the fermion part of the entropy becomes 

T.S a = /* ^ [ - ln(l - exp(-pVif)]. (23) 

^ Jo vr exp{-fiu/K) - 1 



The phonon part S p h is given by the same expression as above but with K = 1. Therefore, 
the entropy S S u(K) of the long-range model is written as S S u(K) — ~d^su(K)l '9T with 

dp 



-T r W ln[l + e-M* + e~ 2 ^l K + ... + e -^-D^]. (24) 

JO 7T 



This is in fact valid for any T. It is natural to interpret ^su(k) as that of ideal para-fermions 
of order K for which up to K — 1 particles can take the same quantum number p. The para- 
fermion is reduced to the fermion in the case of K = 2, which can also be regarded as a pair 



of spinons [17]. In the case of n < 1, the ho Ion also contributes to the entropy and Eq.(f24[) 
describes the spin part at T -C T mix . 

In summary we have derived thermodynamics of the supersymmetric t- J model and have 
shown that the system at low T is equivalent to a set of ideal particles obeying fractional 
statistics. The strong spin-charge separation is caused by the absence of interaction among 
them. We note that description of fractional statistics is the same as that of single-component 
systems only at low T. At higher T the multi-component character of fractional statistics 
appears explicitly. Detailed discussion for general T with extensive numerical results will be 
given separately. The authors thank N. Kawakami for useful discussions. 
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Figure 1. The spin susceptibility against temperature for various fillings. The unit of energy 
is t. 
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